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Abstract 

We consider the discrete time unitary dynamics given by a quantum walk on the lattice 
Z rf performed by a quantum particle with internal degree of freedom, called coin state, 
according to the following iterated rule: a unitary update of the coin state takes place, 
followed by a shift on the lattice, conditioned on the coin state of the particle. We study 
the large time behavior of the quantum mechanical probability distribution of the position 
observable in Z rf when the sequence of unitary updates is given by an i.i.d. sequence of 
random matrices. When averaged over the randomness, this distribution is shown to display 
a drift proportional to the time and its centered counterpart is shown to display a diffusive 
behavior with a diffusion matrix we compute. A moderate deviation principle is also proven 
to hold for the averaged distribution and the limit of the suitably rescaled corresponding 
characteristic function is shown to satisfy a diffusion equation. A generalization to unitary 
updates distributed according to a Markov process is also provided. 

An example of i.i.d. random updates for which the analysis of the distribution can be 
performed without averaging is worked out. The distribution also displays a deterministic 
drift proportional to time and its centered counterpart gives rise to a random diffusion 
matrix whose law we compute. A large deviation principle is shown to hold for this example. 
We finally show that, in general, the expectation of the random diffusion matrix equals the 
diffusion matrix of the averaged distribution. 

1 Introduction 

Quantum walks are models of discrete time quantum evolution taking place on a d-dimensional 
lattice. Their implementation as unitary discrete dynamical systems on a Hilbert space is 
typically the following. A quantum particle with internal degree of freedom moves on an 
infinite d-dimensional lattice according to the following rule. The one-step motion consists 
in an update of the internal degree of freedom by means of a unitary transform in the 
relevant part of the Hilbert space followed by a finite range shift on the lattice, conditioned 
on the internal degree of freedom of the particle. Quantum walks constructed this way 
can be considered as quantum analogs of classical random walks on lattices. Therefore, in 
this context, the space of the internal degree of freedom is called coin space, the degree of 
freedom is the coin state and the unitary operators performing the update are coin matrices. 

Due to the important role played by classical random walks in theoretical computer 
science, quantum walks have enjoyed an increasing popularity in the quantum computing 
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community in the recent years, see for example [28], [3], [21], [23]. Their particular features 
for search algorithm is described in [34], [4], [27] and in the review [31]. In addition, quantum 
walks can be considered as effective dynamics of quantum systems in certain asymptotic 
regimes. See e. g. [11], [1], [28], [26], [9], [30], for a few models of this type, and [6], 
[8], [12], [14], [5] for their mathematical analysis. Moreover, quantum walk dynamics have 
been shown to be an experimental reality for systems of cold atoms trapped in suitably 
monitored optical lattices [19], and ions caught in monitored Paul traps [37]. 

While several variants and generalizations of the quantum dynamics described above are 
possible, we will focus on the case where the underlying lattice is Z d and where the dimension 
of the coin space is 2d. We are interested in the long time behavior of quantum mechanical 
expectation values of observables that are non-trivial on the lattice only, i.e. that do not 
depend on the internal degree of freedom of the quantum walker. Equivalently, this amounts 
to studying a family of random vectors X n on the lattice Z d , indexed by the discrete time 
variable, with probability laws f(X n = k) = Wk(n) defined by the prescriptions of quantum 
mechanics. The initial state of the quantum walker is described by a density matrix. 

The case where the unitary update of the coin variable is performed at each time step 
by means of the same coin matrix is well known. It leads to a ballistic behavior of the 
expectation of the position variable characterized by E^( n )(X n ) ~ nV when n is large, for 
some vector V . This vector and further properties of the motion can be read off the Fourier 
transform of the one step unitary evolution operator. 

In this paper, we consider the situation where the coin matrices used to update the 
coin variable depend on the time step in a random fashion, that is a situation of temporal 
disorder. Let us describe our results informally here, referring the reader to the relevant 
sections for precise statements. 

We assume the sequence of coin matrices consists of random unitary matrices which are 
independent and identically distributed (i.i.d.) and we analyze the large n behavior of the 
corresponding random distribution W u (n) of X%. We do so by studying the characteristic 
function = ¥, W ut n \(e iyX ™). In Section 2, we first show a deterministic result saying 

that the characteristic function at time n can be expressed in terms of a product of n 
matrices, Mj, each Mj depending on the coin operator at step j only, in the spirit of the 
GNS construction, see Propositions 2.9, 2.13. In the random case, the M/s become i.i.d. 
random matrices M w . 

Then we address the behavior of the averaged distribution w(n) = ¥, L0 (W U1 (n)) of X n , 
for n large in Section 3. Theorem 3.10 says under certain natural spectral assumptions on 
the matrices E W (M W ) that X n displays a ballistic behavior 

E w{n) {X n )~nreR d 

where f is a drift vector depending only on the properties of the deterministic shift operation 
following the random update of the coin state. Moreover, the centered random vector 
(X n — nr) is shown to display a diffusive behavior characterized by a diffusion matrix D we 
compute: 

^w(n)((X n ~ nr)i(X n - nr)j) ~ raB^, i, j = 1, 2, ■ ■ ■ ,d. 

We also show in Theorem 3.10 that for any t > 0, y £ M. d , the averaged and rescaled 
characteristic function e~ l ^ ry / ^¥. u (^ t Ay / ^/n)) converges for large n, in a certain sense, 
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to the Fourier transform of superpositions of solutions to a diffusion equation, with diffusion 
matrix H>(v), v G T d , the d-dimensional torus: 

e -.[*/^($^ nl (y/Vn)) -> / e-t^ B ^dv/(27r) d . 

In Section 4, we briefly discuss the relationship between the drift vector r and the 
diffusion matrix D in case the deterministic shift can take arbitrarily large values. Then we 
investigate finer properties of the behavior in n of the averaged distribution in Section 5. 
Theorem 5.2 states that a moderate deviation principle holds for w(n): there exists a rate 
function A* : R d — > [0, oo] such that, for any set r 6 R d and any < a < 1, asn— > cc, 

P(I n -nren( Q+1 »/ 2 r)-e- nam ^ rA ' W . (1.1) 

In Section 6, we consider a distribution of coin operators which allows us to analyze 
the random distribution W u (n), without averaging over the temporal disorder. This dis- 
tribution is supported, essentially, on the unitary permutation matrices. We show that in 
this case W u {n) coincide with the distribution of a Markov chain with finite state space 
whose transition matrix we compute explicitly. Consequently, we get that the centered 
random vector X% — nr converges in distribution to a normal law j\f(0, X), with an explicit 
correlation matrix S, and an explicit deterministic drift vector r given in Theorem 6.6. In 
turn, this allows us to show in Corollary 6.8 the existence of a random diffusion matrix B w 
such that 

E w „ {n) {(X% - nr)i{X» - nr),) ~ nD£, i, j = 1, 2, • • • , d, 

whose matrix elements Wfj are distributed according to the law of X^Xf , where the vector 
X^ is distributed according to Af(0, £). Finally, a large deviation principle for the random 
distribution W w (n) is stated as Theorem 6.14. This example also shows that we cannot 
expect almost sure convergence results for random quantum walks. 

We close the paper by showing how to generalize the results of Sections 3 and 5 to 
the case where the random coin matrices are not independent anymore and are distributed 
according to a Markov process with a finite number of states. See Section 7. 

Let us comment about the literature. In a sense, the situation we address corresponds 
to the cases considered in [29], [17], [15] where the dynamics is generated by a quantum 
Hamiltonian with a time dependent potential generated by a random process. For quantum 
walks, the role of the random time dependent potential is played by the random coin 
operators whereas the role of the deterministic kinetic energy is played by the shift. 

Quantum walks with unitary random coin operators have been tackled in some numerical 
works, see [32], for example. On the analytical side, we can mention [24] (see also [16]) where 
particular hypotheses on the coin matrices reduce the problem to the study of correlated 
random walks. During the completion of the paper, the preprint [2] appeared. It reviews 
and addresses several types of quantum walks, deterministic and random, decoherent and 
unitary. In particular, the averaged dynamics of random quantum walks of the type studied 
in the present paper are tackled, by means of a similar approach. The results we prove, 
however, are more detailed and go beyond those of [2]. 
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We finally note that there exist another instance of random quantum walks in which 
the randomness lies in space rather than in time. In the present context, this means the 
coin operators depend on the sites of the lattice and are chosen according to some law, 
in the same spirit as for the Anderson model. Dynamical or spectral localization are the 
phenomena on interest there. See for example [22], [16], [33] and references therein for results 
about such questions. 

Acknowledgements It is a pleasure to thank L. Bruneau, E. Hamza, M. Merkli and 
C.A. Pillet for fruitful discussions and suggestions about this work. 

2 General Setup 

Let % = C 2d l 2 (Z d ) be the Hilbert space of the quantum walker in Z d with 2d in- 
ternal degrees of freedom. We denote the canonical basis of C 2d by {|r)} re/ d, where 
I± = {±1, ±2, . . . , ±d}, so that the orthogonal projectors on the basis vectors are noted 
P T = |t)(t|, r G I±. We denote the canonical basis of l 2 {'L d ) by {\x)} xl z%d. We shall write 
for a vector ip G H, ip = *}>2 xe zd ip{x)\x) , where i/j(x) = (x\ip) G C 2d and ^2 xeZ d \\ip{x)\\^ 2d = 
||V>|| 2 < oo. We shall abuse notations by using the same symbols (•[•) for scalar products 
and corresponding "bra" and "ket" vectors on U, C 2d and l 2 (Z d ), the context allowing us 
to determine which spaces we are talking about. Also, we will often drop the subscript C 2d 
of the norm. 

A coin matrix acting on the internal degrees of freedom, or coin state, is a unitary matrix 
C G M2d(C) and a jump function is a function r : I± — >■ 7L d . 

The corresponding one step unitary evolution U of the walker on % = C 2d (g) l 2 (7L d ) is 
given by 

U = S(C®I), (2.1) 
where I denotes the identity operator and the shift S is defined on % by 

S = ^ ^2 P T <E)\x + r(T))(x\ + P^ T (g)\x + r(-T))(x\ 
x£Z d re{l,---,d} 

= Pt ® \x + r{r)){x\. (2.2) 

By construction, a walker at site y with internal degree of freedom r represented by the 
vector |t) \y) G % is just sent by S to one of the neighboring sites depending on r 
determined by the jump function r(r) 

S |t) ® \y) = \t) ® \y + r(r)). (2.3) 

The composition by C (g> I reshuffles or updates the coin state so that the pieces of the 
wave function corresponding to different internal states are shifted to different directions, 
depending on the internal state. We can write 

U=Y^ ^2 PrC ®\x + r(T))(x\. (2.4) 
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Given a set of n > unitary coin matrices C k G M2<z(C), = 1, ■ ■ • , n, we define the 
corresponding discrete evolution from time zero to time n by 

U(n,0) = U n U n - 1 ---Ui, where U k = S (C fc ®I). (2.5) 

Let / : Z d — >■ C and define the multiplication operator F : -D(F) — >■ % on its do- 
main D(F) C % by (Fip)(x) = f(x)ip(x), Vx G Z d , where V G -D(F) is equivalent to 
^2 xG id \f (x)\ 2 \\ip(x)\\ 2 ~, d < oo. Note that F acts trivially on the coin state. 

When / is real valued, F is self-adjoint and will be called a lattice observable. 

2.1 Vector states 

In particular, consider a walker characterized at time zero by the normalized vector tpQ = 
tpo <8> |0), i.e. which sits on site with coin state ipo. The quantum mechanical expectation 
value of a lattice observable F at time n is given by (F)^ (n) = (^o\U(n,0)*FU(n,0)^o)- 
As in [16], a straightforward computation yields 

Lemma 2.1 With the notations above, 

U(n,0) = Y, E Pr n C n Pr n . 1 C n - 1 ---P Tl C 1 

x& d Ti,T 2 ,...,T n e/± n 

®\x + r(ri) H h r(r n ))(x| 

= E E «/fc(n)0|x + fc)(x|, (2.6) 

where 

J k {n)= Yl Pr n C n Pr n . 1 C n -l---Pr 1 C 1 eM 2d (C) (2.7) 

">"i,T2,...,T n ei^ n 

EJ =1 r(r s )=fc 

and Jfc(n) = 0, ifJ2™=i r ( T s) / Moreover, for any lattice observable F, and any normal- 
ized vector tpo = ® |0) ; 

(F)^(n) = (Vo|^*(n,0)FC/(n,0)^ o )= J! /(fc)(^o|^(n)Vfc(n)^o) 

fcez d 

= £/(fc)W fc (n), (2.8) 

where W k {n) = \\Jk{n)^Q\\ 2 c2d satisfy 

Y W k (n) = Y ||^H^|lc 2d = WMn = 1- (2-9) 

kez d kez d 

Remark 2.2 VFe view the non-negative quantities {W k (n)} n ^* as the probability distri- 
butions of a sequence of ' 7L d -valued random variables {X n } n& f$*, with 

Prob(X„ = k) = W k {n) = (i> o \U(n,0)*(I®\k)(k\)U{n,0)i;o) = \\Jk{n) Vo \\ 2 c2d , (2.10) 

in keeping with (2.8). In particular, {F)^ (n) = K Wk ^(f(X n )). We shall use freely both 
notations. 
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Remark 2.3 All sums over k G Z k are finite since Jfc(n) = if maxj^...^ \kj \ > pn, for 
some p > since the jump functions have finite range. 

We are particularly interested in the long time behavior, n ^> 1, of (X 2 )^ (n), the expec- 
tation of the observable X 2 corresponding to the function f(x) = x 2 on Z d with initial 
condition tpo. Or, in other words, in the second moments of the distributions {Wk(n)} ne fq* . 

Let us proceed by expressing the probabilities Wk{n) in terms of the C^'s, k = 1, . . . , n. 
We need to introduce some more notations. Let I n (k) = {n,--- ,r n }, where t\ G I±, 
I = 1, . . . ,n and Yl?=i r ( T i) = In other words, I n (k) denotes the set of paths that link 
the origin to k G Z d in n steps via the jump function r. Let us write (po = J2 T el± a r\ T )- 

Lemma 2.4 

n 

W k {n) = J2 VMCI rOinldro) IJ^-il^ ^-i>- (2-H) 

T 0'{ T 1.'" ,Tn}eIn(k) s = 2 

s.t. t u = t^ 

We approach the problem through the characteristic functions $ n of the probability 
distributions {W.(n)} n< =^* defined by the periodic function 

$ n (y) = E w(n) (e^) = Wkinyy", where y G [0, 2vr) d . (2.12) 



To emphasize the dependence in the initial state, we will sometimes write $5?° and/or 
W£°(n). All periodic functions will be viewed as functions defined on the torus, i.e. 
[0, 27r) d ~ T d . The asymptotic properties of the quantum walk emerge from the analy- 
sis of the limit in an appropriate sense as n — > oo of the characteristic function in the 
diffusive scaling 

lim $> n {y/Vn) (2.13) 

Looking for a relation between Wk{n) and Wk{n + 1), we find that the condition r n = T' n 
is a nuisance we can relax now and deal with later. 

Consider the set of paths G n (K) in J? d from the origin to K = ^ Yia ^ ne 

(extended) jump function defined by 



R:I\^Z 2 \ R( T ;) = (lpr]), (2-14) 



that is paths of the form (Ti, • • • , T„_i,T n ), where T s = G I±, s = 1,2, ... ,n, and 

X)™=i R{T S ) = K. Then note that the generic term in Lemma 2.4 reads 



(r' s -i\C* s t' s )(t s \C s r s _i) = (r^|C s t^X^C, r s _i) = foO^KCoC.) t^iOt^), (2.15) 
where, in the last expression, we introduced the unitary tensor product 

V(s) = C S ®C S in £ 2d ®£ 2d . (2.16) 
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The canonical basis of C 2d ® C 2d is {|r t')} ^ TT ,^ eI 2 , hence, with the identification 

T=(1\~\t®t'), (2.17) 
we can write the matrix elements of V{s) as 

(<j®a'\(C s ®C s )T®T')^V(s) S T, S=fy,T=fyell. (2.18) 
With the decompositions 

= ^ a T \r) Xo = <£o®W = Y a ra T '\r <g> r) (2.19) 

t£/± (t,t')£/| 



we can write 



(HlF^xo} = V{l) TlTo A To , where A To = a To a T ,, for T = Q 



(2.20) 



With these notations, we consider the weight of n-step paths in Z 2d from the origin to 
K, with last step T, defined by 

W T K {n)= Yl V(n) TTn _ 1 ---V(2) T2Tl V(l) TlX0 . (2.21) 

(Tl,-'-,T n _ 1 )el2 n ^ 1 s.t. 
.T , n _ 1 ,T)SG„(X) 

Note that by construction, see Lemma 2.4, 

W fc (n) = ^ W^ fc (n) with = Q and ff ± = ,r € J ± } . (2.22) 

We also introduce corresponding periodic functions n > by 

$£(y) = e iYX w£(n), where Y G T M (2.23) 



and, see (2.20), 

^(y) = A T . (2.24) 
These definitions lead to the sought for relationships: 
Proposition 2.5 For all n e N, T e 1±, K G 1 2d and Y G T 2d , 

W£(n + 1) = Y. V ^ n + 1 >sW S K _ R{T) {n) : (2.25) 
Sell 

< + i(Y) = Y et¥RiT)v ( n + l ^s<(Y). (2.26) 



We can express these relationships in a yet more concise way as follows. Recall that for 
each basis vector T ~ r ® r' and each vector Y = (y, y') € T d x T d , we have 

YR(T) = yr(r) + yV(r') G R. (2.27) 

Introduce the vectors in C 4d2 ~ C M <g> C M with y G T M and n > 

*nP0= ^ ^(yJlT®/) and * = ^tIt^t') (2.28) 

T=(r,r')e4 T=(r,r')G4 

where <&o is determined by the internal state (po only and is independent of Y. 
With the matrices on C 2d <8> C 2d ~ C 4d expressed in the ordered basis I± by 

D{Y)= k®T y )(r®r'|, with y G T M , (2.29) 

T=(r,r')64 

and 

^) = c s ®c si M s (y) = D(y)y( s ), (2.30) 

we get the 

Corollary 2.6 For any n > and Y G T 2d , 

* n (y) = M n (y)M n _i(y)---Mi(y)* - (2.31) 

Remark 2.7 TTie matrix D(Y) can be expressed as a tensor product of unitary diagonal 
matrices. Let Y = (y,y') G T d x T d . Then 

D(y,y') = Yj e i{yr{T)+v ' T{T,)) \T®T'){T®T'\ = ^ e ij/r(r) |t>(t| <8> e^'VX^I 

(r,r')e/| rt/± r'e/± 

= d(y)®d(y') where d(y) = ^ e ij/r(r) |r)(T|. (2.32) 

r6i± 

Consequently, we can write 



M s (y,y') = d{y)C s ® d(-y')C s . (2.33) 
Together with the fact that <&o = <A) ® Wo> this yields 



* n (y, y) = d(y)C n • • • d(y)Ci y> ® d(-y')C n d{-y') ■■■C 1 ip 

= Jn(y)vo ® J n (-y')<po = (Jn(y) ® J n {-y')) <P0 ® W- ( 2 - 34 ) 

We note here for future reference that i7 n (y) is the Fourier transform of Jfc(n): 
Lemma 2.8 For any y 6 T d 

J" n (y) = ^ e^ fe J fc (n). (2.35) 



Proof: With the convention (2.27), the right hand side reads 

fcez d (r 1 ,T 2 ,...,T„)e/ ± " 

^ e ^(m)| Tn)(rn | Cne ^(r„_ 1 )| rn _ l)(Tn _ 1 | C . n _ 1 . . . e^lnXnld 

(ri,T 2 ,...,T n )e/ ± n 

= d(j/)C„d(j/)C n _i • • • d(j/)Ci = J n (y). (2.36) 



Eventually, the characteristic function <3?„(y) we are interested in, see (2.12), can be 
obtained from § n (Y). We shall denote the normalized measure on the torus T d by dv = 

dv 

Proposition 2.9 For any n G N and y G T d , 

$E°(2/)= / <*i|M n (y-t;,i;)M n _i(j/ -«,«)■■ -Mi (j/-t;,t;)*o)dw, (2-37) 

w/iere 

*i= E i T >= E i r ® r >- ( 2 - 38 ) 

TeH ± r6/± 

Proof: Following (2.22), to get <& n (y) from {$^(y)} TgJ 2 we need to restrict the sum in 

(2.23) to K = K k , k G Z d , and to sum on all last steps T G H±. 
Let a be the distribution on C°°(T d x T d ) defined by 



«(/(■,■))= / f(v,-v)dv. (2.39) 
Its Fourier coefficients satisfy ) = for all if = (&, fe') G l? d so that 

2 e^W&n)^* = a * *£(y) = a(#£(y - •)). (2.40) 

K=(k,k') 

On concludes using periodicity and by observing that the form yields the summation 
on T G H±. ■ 

Remark 2.10 Noting that *i = J2 T ei± \ T ® r )> we S et 

= [ Tr{J n {y-v)\w){w\Jn(.-v))dv 

= [ Tr(JZ(-v)J n (y-v)\<po)(<Po\)dv. (2.41) 
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2.2 Density matrices 

The analysis above can easily be adapted in order to accommodate more general initial 
vectors or density matrices. 

A density matrix p is a trace class non- negative operator on % = C 2d <X> l 2 (7L d ) which 
can be represented by its kernel 

p = (p(.x,y)) {x y )eZdxZd d, where p(x,y) G M 2( ^(C) (2.42) 

such that 

p= P{^y)®\x){y\. (2.43) 

The matrix p(x, y) satisfies 

p(x,y) = p*(y,x) => p(x,x) = p*{x,x) > (2.44) 
and its elements are given by 

Pa,r{x,y), (a,T)eI±, so that (a <g> x\p t <g> y) = p a , T {x, y) ■ (2.45) 
Since p > is trace class, and C d is finite dimensional, we have 

Trp(x,x) = ||p||i < oo, (2.46) 

\\p{x,x)\\ < Tvp(x,x) < 2d||p(ar,x)||. (2.47) 

The expectation value of a lattice observable F = 1 f in the state corresponding to p 
reads 

(F) p = Tr(p(I ® /)) = Y, /(^)Tr(p(x, x)), (2.48) 

where the first trace is on % and the second on C 2d , assuming that the sum converges. 

If po denotes the initial density matrix, its evolution at time n under U (n, 0) defined by 
(2.5) is given by 

p n = U{n,0)p o U*(n,Q) (2.49) 
and the expectation of the lattice observable F is denoted by 

(F) PQ (n) = Tv( Pn (l0f)), (2.50) 

if it exists. 

Let us specify regularity properties on the lattice observable F = I (g) / and the initial 
density matrix po which imply that all manipulations below are legitimate. 

Assumption R: 

a) The lattice observable is such that, for any p, < oo, 3C P < oo such that 

\f(x + y)\<C„\f(x)\, V(x,y) € Z d x Z d with ||y|| < p. (2.51) 
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b) The kernel po(x, y) is such that 

J2 \\po{x,y)\\ <oo (2.52) 

(x,y)eZ d xZ d 

\f(x)\\\po(x,x)\\ <oo. (2.53) 

xez d 

In a similar fashion to Lemma 2.1, we can express (F) po (n) in the following way 
Lemma 2.11 The kernel of p n reads 

p n (x,y)= ^2 J k (n)p (x - k,y - k')J£,(n). (2.54) 

0,fc')ez d xz d 

Let F = I (g) / and />o satisfy Assumption R. Then p n satisfies Assumption R b) and 
(F) P0 (n) = E MJk(n)po(z-k,z-k')J* k ,(n)) 

zeZ d (k,k')&Z d xZ d 

= E^) E 1t(^(n)J fc (n)po(z-A;,z-fc')). (2.55) 

zez d (k,k')ez d xz d 

Proof: Since for all n G N, the summations on k and /c' are restricted to ||A;|| < p(n) and 
||</fc( n )|| < c(n), we need to control 

S : = El/^l E \\Po(z-k,z-k')\\ (2.56) 

zez d (k,k')ez d xz d 

under Assumption R. Now, po > implies P xy PoPxy > 0, where P xy is the orthogonal 
projector on % 

P xy = I®\x)(x\ + I®\x)(y\ + I®\y)(x\ + I®\y)(y\. (2.57) 
In other words, the following Ad x Ad block matrix is non-negative 

p(x,x) p(x,yf 



p*(x,y) p(y,y)J' 

According to Lemma 1.21 in [36], this is equivalent to p(x,x) > 0, p(y,y) > and 3W, 
\\W\\ < 1 such that p(x,y) = p(x,x) l / 2 W p{y,y) 1 / 2 . Hence, 

\\ P {x,y)\\<\\p{x,x)tl^ P {y,y)tl 2 . (2.59) 

Applied to (2.56), this yields together with Assumption R and Cauchy Schwarz, 

S < E I/WI 1/2 I/WI V2 E \\po{z-k,z-k)\\V 2 \\p Q {z-k',z-k')\\V 2 

zez d (k,k')ez d xz d 

< c »(n) E ^(mz-mipoiz-Kz-m) 1 / 2 



X 



(fe,fc')ez d xz d 

||fc||<M("),l|fc'll<M(«) 



X 



(\f(z-k>)\\\p (z-k>,z-k>)\\) 1 / 2 



< C Kn) \ Y, X | J2\f(x)\\\po(x,x)\\<^. (2.60) 

||ifc||</i(n),||fe'||</i(n) 
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This generalization of Lemma 2.1 allows us to give an interpretation in terms of classical 
random walk on Z rf 

P(X n = z) = Tr(p n (z,z)) = W z Po (n), (2.61) 

with corresponding characteristic function <&n{y) = Y^z&z d e tyz Wz°{n). 

Due to the expression of W z (n) as a convolution, the characteristic function will be 
expressed product of Fourier transforms. 

For Y = (y, y') G T d x T d , we define the matrix valued Fourier transform of a density 
matrix pQ by 

Ro(Y)= Yl e^ k+ y'^p (k,k')eM 2d (C). (2.62) 

(k,k')eZ d xZ d 

Because of (2.46) and (2.44), Rq is uniformly continuous in Y and satisfies 

R* (y,y') = Ro(-y',-y). (2.63) 

Then, the Fourier transform of J n being J Lemma 2.8, the Fourier transform R n of p n 
reads 

Rn(y,y') = £ e^+^') Y Jk(n)po(x-k,x'-k')J* k ,(n) 

(x,x')eZ d xZ d (k,k')eZ d xZ d 

= My)Ro(y,y')J:(-y') (2-64) 

Proceeding as above, we arrive at the generalization of (2.41) 
Lemma 2.12 For any y G T d , 

K°(y) = I Tr R n (y-v,v)dv (2.65) 

JT d 

= [ Tr(J n (y-v)Ro(y-v,v)J*(-v))dv 
JT d 

Tr (J*(- v )J n {y - v)R (y -v,v))dv. 



/ T d 

Let 

Ro(2/,z/)= £ (r\Ro(y,y')r')\r^r'} eC M «C M . (2.66) 

(r,r')€4 

Then, making use of the identity 



M n (y, y')M n _i(y, y') ■ ■ ■ M^y^y 1 ) = J n (y) J n (-y'), (2.67) 
it is straightforward to get from (2.65) the following generalization of Proposition 2.9 
Proposition 2.13 

K°(y)= [ (*i\M n (y-v,v)M n - 1 (y-v,v)---M 1 (y-v,v)R (y-v,v))dv. (2.68) 

JT d 
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Remark 2.14 The map y i— >■ Ro(y — i>, u) is continuous only under (2.46). Under As- 
sumption R for an observable increasing at infinity, this map becomes more regular. 

Remark 2.15 The procedure consisting in extending the space C 2d where the Cj 's act to 
the tensor product C 2d <S> C 2d where we consider Cj Cj is parallel to the GNS construction. 
Once the Fourier transform in position space is taken, it allows us to write the action of 
the one-step dynamics in the coin variables as the action of a unitary matrix in C 2d C 2d 
and to replace the density matrix by a vector. 

3 Random framework 

For a deterministic non periodic set of coin operators, not much can be said about (F)^ (n) 
in general. Therefore we consider the following random quantum dynamical system which 
defines a quantum walk with random update of the internal degrees of freedom at each 
time step. Let C(lo) be a random unitary matrix on C 2d with probability space (f2,cr, d/x), 
where d/j, is a probability measure. We consider the random evolution operator obtained 
from sequences of i.i.d. coin matrices on (S7 N * , J 7 , cflP), where T is the a-algebra generated 
by cylinders and dF = <8>fc g N*d/x, by 

U w {n, 0) = C/ n (uJ)£/ n _i(uJ) • • • C/i(aJ), where U k {oJ) = S (C(oj k ) ® I), (3.1) 

and uJ = ((jJi,(jJ2,0J3, . . . ) G . The evolution operator at time n is now given by a product 
of i.i.d. unitary operators on %. We shall denote statistical expectation values with respect 
to P by E. 

All results of the previous section apply, with each occurrence of C s replaced by C(oj s ). A 
superscript uJ will mention the resulting randomness of the different quantities encountered. 
In particular, the random dynamical system at hands yields random matrices J k ( n ) £ 
A^2cz(C), which, in turn, define random probability distributions {W£(n)} n€ ^* on Z d which 
satisfy (2.9) for all n G N* and oJ G The corresponding characteristic functions ^ 

become random Fourier series whereas &n(Y) 1S obtained by the following product of i.i.d. 
random matrices 

*£(y) = M Wn {Y)M u)n _ 1 (¥)■■■ M m (r)* (3.2) 

where, with Y = (y,y'), 

(Y) = D(Y)V(u s ) = d(y)C(u s ) ® d(-y')C(uj s ) (3.3) 
are distributed according to the image of dfi by the inverse mapping C i-» D(Y)C C. 

3.1 Diffusive averaged dynamics 

We consider in that section the statistical average of the motion performed by random 
quantum walk, and, more specifically, its diffusive characteristics. For the lattice observable 
X 2 , we will derive results regarding the long time behavior of 

E((A 2 )^ o )(n) = E([/±(n,0)Vo|^ 2 ^(n,0)^ > = E(E^ r(n) (X T 2 J). (3.4) 
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It means, see (2.8), that we consider the motion corresponding to the averaged probability 
distributions defined by 

w k (n) := E(Wf (n)), k G Z d , n G N*, (3.5) 
with corresponding characteristic function 

<Mj/) = E w(n) (e^) = ^ «, fc (n)e ifc ». (3.6) 

Remark 3.1 To stress the dependence on uJ in the distribution W£(n), we shall denote 
the corresponding random vector on the lattice by Xf L . When we consider the averaged 
distribution w(n) instead, we shall write X n for the corresponding random vector. 

Our results generalize those of [24] in the sense that the distribution of the random coin 
matrices considered here is arbitrary. As a consequence, the analysis cannot be mapped to 
that of a persistent or correlated classical random walk on the lattice, as was observed in 
[24]. 

Let £ and M(Y) be the matrices defined by 

£ = E(V(uj)) =E(C(uj)®C(uj)) and M{Y) = D(Y)£. (3.7) 

Note that while V(lu) is a unitary tensor product, its expectation £ is neither unitary, nor 
a tensor product in general. But ||£|| < 1 and therefore ||.M(y)|| < 1. 

Since the {C(w s )} s6 n* are i.i.d., we immediately get from Corollary 2.6 and Proposition 
2.9 that 

E(*£)(y) = (M(y)) n * , (3-8) 

so that 

HK){y)= I (*i\(M(y-v,v)) n * )dv. (3.9) 

The analysis of the diffusive scaling limit (2.13) now relies on the spectral properties of the 
matrices £ and M(Y). 

3.2 Spectral properties 

The structure of these matrices implies the following deterministic and averaged statements: 

Lemma 3.2 Let V(u) = C(u)0C(w), £ = E(V(w)), M(Y) = D{Y)£ and let S denotes 
the unitary involution defined by Sip tp = tp ® ip, for all tp, tp € C 2d . 

Then, for all oj and all y, *i = Yl T ei± \ T ® T ) ^ s invariant under V(uj),£,Ai(y, —y),S 
and their adjoints. Consequently 

\\£\\ = Spr(£) = \\M(y,-y)\\ = Spr (M(y, -y)) = 1, (3.10) 

where Spr denotes the spectral radius. Moreover, 

SV(co)S = V(oj), SM(y,-y)S = M(y,-y), S£S = £, (3.11) 

so that 

a{V{u)) = WW), <£) = W), <r(M(y, -y)) = a(M(y,-y)). (3.12) 
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Remark 3.3 For MiY), we only have 

Spr(M(y))< ||M(y)|| < 1. (3.13) 

Remark 3.4 Before taking expectation values, 1 is a eigenvalue at least 2d times degenerate 
for the unitary matrices V(oj) and M U! (y, —y), for all w and y, because of their tensor product 
structure (2.16), (2.33). 

We shall work under an assumption which implies that in the long run, the averaged 
quantum walk loses track of interferences and acquires a universal diffusive behavior. At 
the spectral level, this is expressed by the fact that after taking expectation values, 1 is the 
only eigenvalue of M(y, —y) = E(M w (y, — y)) on the unit circle and it is simple. 

Let D(z, r) C C denote the open disc of radius r centered at z G C. 

Assumption S: For all v G [0, 2ir) d = T d , 

a(M(—v,v)) fl dD(0, 1) = {1} and the eigenvalue 1 is simple. (3-14) 

Remark 3.5 Actually, because of the form of (3.9), it is enough to consider the spectrum 
of the restriction of Ai(—v,v) to the M.*(Y)-cyclic subspace for generated by Set 

1 = Span {AT (y) fc *i, k G N, Y G T d x T d }, (3.15) 

and let Pi = P£ be the orthogonal projector onto 1. If Pi / I, we can work under the 
weaker 

Assumption S': For all v G [0, 2ir) d = T d , 

cr(A4(— v, v)\i) n dD(0, 1) = {1} and the eigenvalue 1 is simple. (3.16) 

Indeed, note that M*(Y)P X = P X M* {Y)P X so that, at the level of linear forms 

(*i|M(r) n = (Ar(y) n Px*i| = (^ 1 \(p x M(Y)p x ) n = <*i|M(y)| x n . (3.17) 

While it is often necessary in applications to use S', see the examples, we keep working 
under S below in order not to burden the notation. 

Let \I/o = ^i/H^i ||. Under assumption S, V&o spans the one dimensional spectral 
subspace of A4(— v,v) associated with the eigenvalue 1. Moreover, by Lemma 3.2, the 
corresponding rank one spectral projector reads P = | 1 I r o)(*o| and is w-independent. With 
Q = I — P, we have the spectral decomposition 

M(-v,v) =P + QM(-v,v)Q (3.18) 

where, under assumption S, 3 e < 1, independent of v G T d such that Spr QA4(— v , v)Q < e. 

In keeping with (3.9) and the diffusive scaling (2.13) to be used below, we perform a 
perturbative analysis of the spectrum of M(y — v,v) for small values of ||y||, uniformly in 
v G T d . Let us introduce the following notation for (y, v) G T d x T d 

M v {y) = M(y - v,v), so that M v (0) = M(-v,v). (3.19) 
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Now, with D(y,y') = d{y) d(y'), see (2.32) 

M v (y) = D(y, 0)M V (0) = M v (0) + ^ (e iyr ^ - l)|r)(r| 8lM„(0) 

rG/± 

= M,(0) + F(j/)M,(0), (3.20) 

where ||^W^ (0) || = 1 and ||.F(y).M„(0)|| < c\\y\\, with c independent of v. 

Since the map (y,v) i-> .M^y) is actually analytic in C ' x C d , we can say more. For 
v > 0, let 7? = {2 G C d | ffiz G T d , Qfy- < 1/, j = 1, . . . , d.} C C d be a complex neighborhood 
of T d . For y > 0, let 73(0, y ) = {y G C d \ \\y\\ < y }- 

Analytic perturbation theory, see [18], then yields the following 

Lemma 3.6 Under assumption S, there exists < S < 1, v = u{5) > and yo = yo(5) > 
such (y,v) £ (7^nB(0,?/o)) x T d implies 

a(M v (y))nD(l,6) = {X^v)} (3.21) 
<r(M v (y))\{Xi(y,v)} C D(0,1-S). (3.22) 

Moreover, \i(y, v) is simple, analytic in (7^ d n£>(0, yo)) x T„ and Ai(0, v) = 1 /or a// v 6 7^ d . 
TTie corresponding spectral decomposition reads 

M v {y) = X 1 (y,v)P(y,v)+M Q (y,v), (3.23) 

w/iere P(y, v) is analytic in (T„ n £>(0, yo)) x and P(0, v) = P for all v G 7^ d . 

W^i/t Q(y,v) = I — P(y,v), the restriction Aiq(y,v) = Q{y,v)A4 v (y)Q(y,v) satisfies 

Spr (M Q (y -v,v)) <l-5. 

We need to compute Ai(y,v) = Tr (P(y, v)M v (y)) to second order in y. We expand 
F(y) as 



F(y) = F 1 (y) + F 2 (y) + 0(\\yf) (3.24) 

1 3 \ 



rel± rel± 

and introduce the (unperturbed) reduced resolvent S v (z) for v G and z in a neighborhood 
of 1 such that 

(M v (0) -z)- 1 = J^ + S v (z) with P= |*o) (*o|- (3.25) 
1 — z 

We have for a simple eigenvalue, (see [18] p. 69) 

Ai(y,t;) = l + TriFiWMvMP) (3.26) 
+ Tr(F 2 (y)M v (0)P - Fi(y).M„(0)^(l)F 1 (y).M„(0)P) + O v (\\yf). 

Explicit computations with symmetry considerations yield 
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Lemma 3.7 For all v G T d and y G B(0,yo), there exists a symmetric matrix B(d) G 
M^(C) suc/i that 

Ai(y,«) = 1 + ^ E^ + ^dl^ll 3 ) 

+ ^(E^Y^+ E (F(r))(F(r'))((T®r|^(l)r'0r')-i} 

= i + ^E^-^^^+^diyii 3 )- ( 3 - 27 ) 

T6/± 

TTie mop f I—?- is analytic in T d ; when v 6 T rf ; D(v) G Md(IR) is non-negative and 

°( v )j,k = %~Wr^(°'' u )> 3,k G {1,2, ... ,d}. Moreover, C^dlyH 3 ) is uniform in v £ T d . 

Proof: Existence and analyticity in v of D(v) follow from analyticity of Ai in y and 
analyticity of S v (l) in v, see (3.25). Since H)(v)j^ = Qy Q yk ^(0; v )i the matrix is symmetric. 
For v G T d , the symmetry (3.11) implies SS V (1)S = S v (l) so that 



(t' ® t'\S v (1)t ® t) = (t' r'|5^(l)5r ® t) = (r' ® t^S^t r). (3.28) 

Hence the matrix elements D(v) for -u G T d are real as well. Finally, (3.13) implies that 
(y\B(v)y) > for all y G T d . ■ 

Remark 3.8 Using the notation f = ^ X^re/± /( r ) ^ or an 2/ function on I±, B>(v) reads 

D(u) = 2|r)(r| - J |r)(r| - I |r(r))(r ® t|^(1)t' ® r')^')!, (3.29) 

z a 

T,r'ei± 

where the "bra" and "ket" notation is understood in R d for vectors r(r) and in C 2d <8> C 2d 
for t r. 

We are now set to prove the 

Proposition 3.9 Under assumption S, uniformly in v G T d , in y in compact sets of C 
and in t in compact sets ofW^, 

lim M^My/n) = e ity¥ P, (3.30) 

n— >oo 

lim MMfo/v^e-*^/^ = e -i {ymv)y) P. (3.31) 

n— >oo 

Proof: Let v G 7^, i G G C M+, and y € K C C, G and if compact. We consider n large 
enough so that y/>/n and y/n belong to B(0,yo), uniformly in y G K. The decomposition 
(3.23) implies for any n large enough, 

Mt\y) = X [ ! n] (y,v)P(y,v)+Ml n \y,v) (3.32) 
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where Lemma 3.6 implies the existence of c > and 1 > 6' > 8, uniform in (y, v, t) £ 
T„ x K x G, such that 

\\M% ] (y,v)\\<c6> tn . (3.33) 



Moreover, by Lemma 3.7, 

(— \ [* n l 

1 + i !l + O v (\\yf/n 2 )j " (3.34) 



n— >oo „ityr 



[tn] 



K (y/Vn,v)e ^ = | e ^M +l _ + o /__))}> (3.35) 



and both P{y/^/n,v) and P(y/n,v) tend to F as n -)> 00, uniformly in (f,y) G 7^ x K. I 

With these technical results behind us, we come to the main results of this section which 
are the existence of a diffusion matrix and central limit type behaviors. 

Let A/"(0, £) denote the centered normal law in R d with positive definite covariance 
matrix £ and let us write X w ~ AA(0, S) a random vector 6 M rf with distribution 
AA(0, X). The superscript uj can be thought of as a vector in R d such that for any Borel set 
AcR d 

POT € A) = ; V I e-^^du. (3.36) 

The corresponding characteristic function is &^(y) = K(e iyXU ') = e~2^l s ^). 

The first result concerning the asymptotics of the random variable X n reads as follows. 

Theorem 3.10 Under Assumption S, uniformly in y in compact sets of C and in t in 
compact sets of~R*^, 

lim $^ (y/n) = e^ (3.37) 

lime"' W ^$r,(^)=/ e -5<»l D ( w )f> d«, (3.38) 

where the right hand side admits an analytic continuation in (t,y) G C x C 2 . 
In particular, for any (i, j) G {1, 2, . . . , d} 2 , 

lim (X ^ o(n) = n (3.39) 

n— >oo fi 



lim ((X-nr)i{X-nr)j)^(n) 



/ Bij^cZv. (3.40) 



n— >oo rj 

Remark 3.11 We ure// caZZ diffusion matrices both B)(v) and D = J* Td B(t>) eft). 
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For any s = (si, s 2 , . . . , s d ) G N d ratt |s| = £? =1 «j and = 1 ■ ■ ■ (j^) d , 

(XfX 2 S2 ---X^)^,(n) _ S2 
Um ((X-nr)f(X- nr)* 2 • • • (X - nr)*<% (n) 



n— »oo n 



1*1/2 



= H) N / (i^e-^W^))!^^, (3.42) 



which shows that all odd moments (\s\ odd) of the centered variable are zero whereas all 
even moments can be computed explicitly. 

Proof: This is a direct consequence of Proposition 3.9 and definition (3.9). The unifor- 
mity of the convergence in the variables (v,y,t) in compact sets provides analyticity after 
integration in v G T d and commutation of the limit and derivations. H 
For initial conditions corresponding to a density matrix po, we have 



Corollary 3.12 Under Assumption S, for any t > 0, 

■ [tn 



hm $f»(y/ n ) = e^, (3.43) 



where 



lim e~* n] vi$f° (y/V^) = / e-3<wl D (")»> (tfilRot-w.v))^ 

= /" e -!(fl D (^)lr( J R (-'y,«))^, (3.44) 
Jf d 



R (-v,v)= e ivl p (k,k + l). (3.45) 



Remark 3.13 Under Assumption R for the observable X 2 , we deduce that for any £ 
{1,2,... ,d} 2 , 

lim {Xi) ^ H) = r u (3.46) 

n— >oo 7j 

hm ^-"^-%(") = / B^Tr^oH^))^. (3.47) 

n->oc n Jjd 



From Corollary 3.12, and Theorem 3.10, we gather that the characteristic function of the 
centered variable X n — nr in the diffusive scaling T = nt, Y = y/y/n, where n — > oo, 
converges to 

/ J / w W2 /XTW t (V)^m-v,v))dv, (3.48) 

where the function under the Fourier transform symbol J 7 is a solution to the diffusion 
equation 

i,j=i x * x i 
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As explained in [17], [15], it follows that the position space density Wk([nt])5{yfnx — k) con- 
verges in the sense of distributions to a superposition of solutions to the diffusion equations 
(3.49) as n — > oo. 

In case where the diffusion matrix B(-u) = D is independent of v, Theorem 3.10 with 
t = 1 says that the characteristic function of the rescaled variable (X n — nr) / yjn defined 
by E u , fe(n )(e jy ( Xn_nr )/ v/ ™) converges to which is the characteristic function of the 

normal law AA(0,D). Hence, by Levy's continuity theorem, see e.g. Theorem 7.6 in [7], 

Corollary 3.14 Assume S, suppose D(t> ) = D > is independent of v G T d . Then, for 
any initial vector = <A) ® 1 0) > we have as n — > oo, in distribution, 

X w ~ A/"(0,B). (3.50) 

Remark 3.15 All results of this section hold under Assumption S' only, mutatis mutandis. 
In particular, if the invariant subspace X coincides with span {\a ® cr)} CT e/ ± ; the matrix 
A4 v (y) is actually independent of v, because D(—v,v) acts like the identity on the latter 
space. Consequently, the diffusion matrix is independent of v as well. 

Remark 3.16 We have chosen to randomize the coin state updates only, but it is possible 
to adapt the method to deal with random jump functions as well. 



Jn 



3.3 Example 



For d = 1, consider the set of three unitary matrices in C 2 given by < I, 



'0 1\ 1 

i oj'Vzyi -i 

and the distribution which assigns the probability p/2 > to the first and second matrices 
and q = 1 — p to the third one. Let r be the jump function defined by r(±l) = ±1 so that 
r = 0. Then, in the ordered basis {| — 1<8>— 1),|1<8>1),| — 1<8>1),|1<8>— 1)}, the corresponding 
matrix £ reads 

(\ 1 q q \ 

1 1 -q -q 
q -q p-q 1 
\q -q 1 p-q) 



(3.51) 



and D(—v,v) = diag(l, l,e , e ). We introduce the following orthonormal basis whose 
first vector is *o : 

{(|-1®-1> + |1®1»/V2,(| 



1 ® -1) - |1 ® 1))/V2, | - 1 ® 1), |1 ® -1)} 



In this basis, M v (0) = D(—v,v)£ writes 

/l 






(3.52) 



M v (0) 






J2v q 







V2 
-ilv q 
V2 





q 

,i2v (p-q) 

2 

„-i2v 1 





q 

V2 
J2vl 
e 2 

-i2v (p-q) 

2 



= i e n v , 



(3.53) 
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where N v is the restriction of M. v (0) to the subspace orthogonal to C\Ev in order to make 
computations easier, we specialize to the case p = l/\/2, so that 

q/V2 = ( P -q)/2=(V2-l)/2 = j. (3.54) 

This way we can write in (C^o)" 1 

/ -1 7 7 \ 

N v - 1 = e i2v 1 e i2v j - 1 e i2v \ . (3.55) 
y e -i2« 7 g -i2«i e ~i2v 1 _ l j 

We have det(iV„ - I) = 2cos(2w)(7 2 + 7) - (2 7 3 + 3/4) < for all v G T so that C* is 
the only invariant subspace under M v (0). Hence S v (l) = (N v — I) -1 . To get the diffusion 
constant I$(v) we need to compute (r (X 1 t|S'(1)t / (8> r') for r, r' = ±1, where S(l) is defined 
on (C^o)" 1 = QC 4 . We have Q| ± 1 ® ±1) = so that 

(t ® t|5(1)t' ® t') _ £ - ,,-VO - ^ ^fg^ . (3-) 

Taking into account the r(r) = r in formula (3.27), we get 

B(v)y 2 = -y 2 (l + 4(^\(N V - l)' 1 ^)), (3.57) 

which, plugging in the value of 7, eventually yields 

= / l 6 - 9 ^+2eo s (2„)(5-4y/2)\ 
V 5V2-4-2cos(2v) / 



4 Einstein's Relation 

An interesting feature of the previous results is that the asymptotic averaged velocity r 
depends on the jump function r only and is independent of the coin distribution. This is 
reminiscent of the asymptotic velocity v(F) reached by a particle subject to a deterministic 
force of amplitude F in a random dissipative environment modeled by random forces. Given 
an asymptotic velocity, the mobility vector fi is defined as the ratio 

u = lim v F /F. (4.1) 

F-5>0 

This mobility (i is then related to the fluctuations of the system around the asymptotic 
trajectory by Einstein's relation which says that the diffusion matrix is proportional to 

IMI- 

In the present framework, neither dissipation nor forces can be directly traced back to 
describe the asymptotic motion 

(X)iB (n) = nr + o(n), n ->■ 00. (4.2) 

Moreover, the motion taking place on a lattice, the asymptotic velocity r has a minimal 
amplitude l/(2d), if it is non zero, which prevents a behavior similar to (4.1). Nevertheless, 
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the jump function r which characterizes the deterministic motion can be thought of as an 
external control parameter, similar to a driving force. 

In order to get an asymptotic velocity which vanishes with the exterior control param- 
eter, we rescale the lattice 2, d to (Z/l) d , with I > 0. This means we consider the variable 

Y n = X n /l e (Z/l) d . (4.3) 

Then we introduce a parameter s £ N as follows. Let r\ and tq be two non-zero jump 
functions such that 

FT = and / 0. (4.4) 
We define a new s-dependent jump function by 

r s (r) = sri(r) + ro(r) £ Z d such that = Yq (4.5) 

and we will consider the large s limit. Hence the rescaled variable Y n satisfies 

lim (1 ^» (n) =^-.= v y (4.6) 

rt-^oo n I 

hm = fV /■ (w) d5 + 0{1/s) ) := Df (4 . 7) 

where ID^ 1 ) (u) is the diffusion matrix computed by means of the jump function ri and the 
remainder term is uniform in v £ T d . Therefore, choosing the scale I = s £ N, we get that 
the diffusion matrix B Y is finite for large s whereas the asymptotic velocity tends to zero. 
Hence, setting F = 1/s, we get for s large 

a = lim v Y /F = r^ (4.8) 
F=l/s->0 



/ B^^dv + oil/s) ^ K\\fi\\, as s oo. (4.9) 



The last formula is admittedly a consequence of a rather ad hoc construction. On the other 
hand, assuming that r± 7^ 0, we get with the same scaling v Y = FT which never vanishes. 

5 Moderate Deviations 

The spectral properties of the matrix M(y — v,v) proven in Section 3.2 allow us to obtain 
further results on the behavior with n of the distribution of the random variable X n defined 
by (2.2). This section is devoted to establishing some moderate deviations results on the 
random variable X n . 

We consider initial conditions of the form V'o = Vo ® 1 0) and we will be concerned with 
X n — nr. Moderate deviations results depend on asymptotic behaviors in different regimes 
of the logarithmic generating function of X n — rir defined for y G M. d by 

A n (y) = \n(E w{n) (e^-™))) G (-00,00]. (5.1) 

This function A n is convex and A n (0) = 0. 



22 



Let {a n } n£ N be a positive valued sequence such that 

lim a n = oo, and lim a n /n = 0. (5.2) 

Define Y n = (X n - nr)/^/na^ and, for any y G R d , let A n (y) = ln(E w ^(e yYn )) be the 
logarithmic generating function of Y n . 

Proposition 5.1 Assume S and further suppose B(t>) > for all v G T d . Let y G R d \ {0} 

and assume the real analytic map T d B v ^ (y\3(v)y) G R+ is either constant or admits a 
finite set {vj(y)}j = i^... t j of non- degenerate maximum points in T d . Then, for any y G M. d , 

lim — A n (a n y) = hyM Vl (y))y) (5.3) 
which is a smooth convex function ofy. 

Proof: This proposition essentially follows from Lemmas 3.6 and 3.7 and the asymptotic 
evaluation of an integral. Let b n = \J a n /n s.t. hin^^oo b n = 0. By construction, A n (a n y) = 
A-n(b n y) where, according to Lemma 3.6 and 3.7, there exists 7 > s.t. for n large enough, 

exp(A n (&„y)) = J (l + f(y\B(v)y) + O v (b 3 n y 3 )) U + O v (b n y)) dv + 0(e^ n ) 

= f e^ {ym ' u)y)+0 ^ anbny3) {I + O v (b n y)) dv + 0{e~^ n ). (5.4) 
Jf d 

All remainder terms O v (---) are analytic in v G 71/, as well as IED(f). An application of 
Laplace's method around each of the non-degenerate maximum points, yields for 1/3 < 
a < 1/2 

exp(A n (6 n y)) = ^e^«^^^^ + °( b ^ 3 ))(G J ( y )/a^ 2 + 0(b n y) + Oa/a^ 1 )) 

+0{e- rn ) + 0{e- Ka ^ 2a ), (5.5) 

where Gj(y) > 0, K > 0, from which the result follows. The case where B(v) is indepen- 
dent of v follows directly from (5.4). The convexity of the limit follows from the convexity 
of A n . The assumed non-degeneracy of the maximum point ensures that the functions 
M d \ {0} B y 1 — y Vj(y) are all smooth by the implicit function theorem. M 

Let us recall a few definitions notations. A rate function I is a lower semicontinuous 
map from M. d to [0, 00] s.t. for all a > 0, the level sets {x | I(x) < a} are closed. When the 
level sets are compact, the rate function / is called good. For any r C l d , T denotes the 
interior of T, while T denotes its closure. 

As a direct consequence of Gartner-Ellis Theorem, see [13] Section 2.3, we get 
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Theorem 5.2 Define A*(x) = sup ygRd ((y\x) - ±(y\B(v 1 (y))y)) , for all x G R d . Then, A* 
is a good rate function and, any positive valued sequence {a n } raG N satisfying (5.2) and all 

r c R d 

- inf A* fx) < liminf — ln(¥((X n - nr) G JrTa^T)) 
xer° n^oo a n 

< limsup — ln(¥((X n - nr) G s/na^T)) < - inf A*(x). (5.6) 

Remark 5.3 As a particular case, when D(u) = D > is constant, we get 

A*(x) = Uxp^x). (5.7) 

Remark 5.4 Specializing the sequence {a n } ra6 N to a power law, i.e. taking a n = n a , we 
can express the content of Theorem 5.2 in an informal way as follows. For < a < 1, 

¥((X n - nr) e n( a+1 )/ 2 r) ~ e""" ini *^ A * W . (5.8) 

For a close to zero, we get results compatible with the Central Limit Theorem and for a 
close to one, we get results compatible with those obtained from a large deviation principle. 



Let us come back to the example in section 3.3. The diffusion coefficient D(f ) given in 
(3.58) admits , as a function of v G T, a single non-degenerate maximum at v = where it 
takes the value B(0) = 2y/2 — 1. Thus we get from the foregoing that a moderate deviation 

2 2 

principle holds for this example, with the good rate function A*(x) = ^ffnT = „,» 



6 Example of diffusive random dynamics 

The results obtained so far can be viewed, essentially, as an adaptation to the quantum 
walk dynamics setup of those proven in [29], [17], [15] for the averaged dynamics and as an 
extension of [24], [16]. 

In this section we consider a specific example of measure d/j, on U(2d), the set of coin 
matrices, for which we can prove convergence results on the associated random quantum 
dynamical system (3.1) for large times, in distribution rather than in average. In particular, 
our example shows that almost sure convergence results cannot be expected in general. 

As noted in Section 3.2, the spectra of V(uj) and M UJ (y,—y) lying on the unit circle 
and admitting 1 as a 2ci-fold eigenvalue prevent us from using the same spectral methods 
as above. For the same reason, the results about products of random contractions in [10] 
do not apply. However, the structure of the example at hand allows for a direct approach 
which, eventually, reduces the analysis to that of a central limit theorem for a Markov chain. 
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6.1 Permutation matrices 

Let &2d be the set of permutations of the 2d elements of I± = {±1, ±2, . . . , ±d}. For 
vr G & 2d and 9 = {9j} jeI± G T 2d , define 

C(7r,9) = E e^MK(r))(r| G f/(2d) so that C CTT (^,0) = e s ° <5 CT)7r(r) . (6.1) 

With A(0) = diag(e ie ) and C(tt) = C(vr,0), we can write 

C(7r,e) = A(e)C(7r), (6.2) 

where C(7r) is a permutation matrix associated with 7r. We recall the following elementary 
properties: For any ir, a G ©2d, 

C(I)=I, C*(ir) = C T (ir) = C(ir- 1 ), C(ir)C(a) = C(ira). (6.3) 

Moreover, Birkhoff-Von Neumann Theorem asserts that the set of doubly stochastic ma- 
trices of order n is the convex hull of the set of permutation matrices of order n whose 
extreme points coincide with the permutation matrices. 

The matrices C(tt, 0) allow for explicit computations of the relevant quantities intro- 
duced in Section 2. It is easy to derive the next 

Lemma 6.1 Let r : I± — > 7L d be a jump function. Given a sequence of n permuta- 
tions 7Ti, 7T2, . . . , 7r n , let (t"i, T2, . . . , T n ) G 7™ ^ e ^ e sequence parametrized by t\ given by 
(Ti,7r 2 (Ti),7r 3 (r2), . . . ,7r n (r n _i)), i.e. suc/i tftai 

Tj = (tTj-tTj-i • • •7r 2 )(ri), j = 2,...,n. (6.4) 

Let 0i, 2 , . . . , ra 6e a set o/ phases, Qj = (j), • • • , n {j)). Then, with the convention 
Cj = C(7Tj, Qj), we get for all k G Z d , 

J k (n) = e i ^+-^)| Tn )( 7ri -i( Tl )|, (6.5) 

TlG/_(- s.t. 

and J k {n) = 0, «/E"=i r ( r j) / 

Consequently, the non-zero probabilities H 7 ^ (n) on Z rf read for any normalized internal state 
vector (po and any density matrix po 

W?(n) = || J k {n)M 2 = E IKVi)N>| 2 , (6.6) 

Tiei-|- s.t. 
E7 =1 t-(^)=* 

W£» = Trp n (M) = E E (^(n)!/^ - j, k - j)^ 1 ^)). 
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Note that the sets of phases ®j, j = 1, . . . , n play no role in the computation of expectation 
values of lattice observables. We set t\ = vri(ro) and note 

<P0 = E a ^.ko) => K^r 1 ^!)!^)! 2 = E KI^tu^)- (6-7) 

roG/± ro£l± 

Hence W^{n) = Z TO ei ± KI 2 *E?=ir(r,-),fc so that for F = 1 ® f and ^ = ^ ® |0)<0| 

= E (»)/(*) = E KIV(X>(t,-)). (6.8) 

-roe/± i=i 



Remark 6.2 /n oi/ier words, given a set of n permutations, there is no more quantum 
randomness in the variable X n , except in the initial state. 

Therefore the characteristic functions take the form 

Corollary 6.3 With Tj = (TTjiTj-i • • • 7Ti)(to), for j = 1, . . . , n, 

n°(v) = E e^-^Vol 2 , (6-9) 
r ei± 

K°(V) = E e iy ^=^) f (r \R (y-v,v)T )dd. (6.10) 

tq£I± 

The dynamical information is contained in the sum S n = Y^j=i r ( T j) which appears in the 
phase. The next section is devoted to its study, in the random version of this model where 
the coin matrices are i.i.d. random variables with values in {C(tt, 0), 7r G &2d, © £ Z M }. 

6.2 Random dynamics 

Assume a random variable C(oj) with values in {C(7r,0) G U(2d), (tt,@) G ©2d x T 2d } is 
defined on a probability space (Q,a,du). The foregoing shows that only the marginal a 
defined on the discrete set {C{n) G U(2d), it G &2d\ or, equivalently on {-7T G &2d}, matters 

m(tt) = /x(a(w) = vr) = = C(tt, 6) | 9 G T m }), tt g & 2d . (6.11) 

We shall use the notation a(co) = uj G ©2d and SI = &2d- The corresponding process is 
denoted by uJ = (uji, uo 2 , W3, ■ ■ ■ ) G il N * and dP = ®k&n*dfi. 

Given ipo G C 2d an initial internal state and a random sequence of permutation matrices 
(oji, . . . ,oj n ), the random variable S n (oJ) = Yl^=i T j(P) e ^ d * s the sum of random variables 
Tj(ZJ), j = 1, . . . , n whose properties are given in the next lemma: 

Lemma 6.4 Let tpo = ^7^ a Tol r o) be the initial condition. The path (to, t±, . . . , T n ) is a 
Markov chain with finite state space I± characterized by the initial probability distribution 



Po(t = cr ) = I a, 



2 



TO 
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and by the stationary transition probabilities 

P(a',a) = Prob(r k (cJ) = a|r k _i(aJ) = a'), fc = 2,3,...,n (6.12) 

given by 

pw,°) = M 71 ")^,^')- ( 6 - 13 ) 

The corresponding transition matrix P = (P(a',a)) G M2d(M + ) is doubly stochastic and 

P = E(C t {oj)). (6.14) 

Remark 6.5 The transition matrix P is unitary iff /j,(tt) = S no>n for some -kq. 

Proof: By Lemma 6.1, Tk(uJ) only depends on {^j}j=i,- ,k an d To and is given by Tfc(uJ) = 
Wkijk-i). Hence 

Prob(r fe (aJ) = a\T k -i(uj) = a') = Prob(a; fc (r fe _i) = cr|r fe _i = a') (6.15) 

= p{{u\u{o') = a}) = ^ ^)S a ,n{a'), 

where we used the independence of the oo^. Finally, (6.1) shows that the right hand side is 
the expectation of C t (uj) w.r.t /U. H 

Considering the diffusive scaling (2.13), we are thus naturally lead to investigate the 
large n behavior of the quantity 

-±=S n (uJ) = -Lf^r^)), (6.16) 

i.e. to a functional central limit theorem for the Markov chain (to, ti, T2, . . . ) with finite 
state space I±, initial probability po and transition matrix P. 

There are simple conditions under which a functional central limit theorem holds for 
Markov chains with finite state space, see e.g. in [7]. Let us recall the few basic notions and 
results associated to Markov chains with finite state space, F, characterized by a transition 
matrix P G M\p\(M+) s.t. J2 T eF P{ a i T ) = 1 that we will need below. 

A transition matrix P is irreducible if, Vcr, r G F, 3n G N* such that P n (cr, r) > 0. A 
probability distribution po, considered as a vector in M) F \, is invariant for the transition 
matrix P if P T po = po- If Po is invariant, then the Markov process is stationary, 

Prob((7i,,ri,---)€B) = Prob((75 fc ,7T fc+ i,---)€ J B), VfceN, BcF N . (6.17) 

If P is irreducible, the invariant distribution po is unique and po(r) > 0, Vr G F. If P is fur- 
thermore doubly stochastic, the invariant distribution uq is uniform uq(t) = 1/\F\, Vr G F. 
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In terms of spectral properties, an irreducible stochastic matrix P admits 1 as a simple 
eigenvalue. If it is furthermore doubly stochastic, the uniform vector uo is invariant under 
both P and P*. 



Hence, if we take as initial distribution the uniform measure uq(t) = \/{2d) Vr £ I±, 
which is invariant for the doubly stochastic transition matrix P = E M (C T (u;)), the Markov 
process is stationary. Moreover, 

rel± 

From Thm 20.1 in [7] and its applications page 177, or [25], we have 

Theorem 6.6 Let (po = Y2 T ei± a r\ T ) and po s.t. po(r) = \a T \ 2 . Assume the transition 
matrix P = E((7 T (o;)) is irreducible. Then, limji^oo ^ X)j=i r ( T j(^)) = ^ almost surely 
and, with convergence in distribution, 



3 

provided the covariance matrix 



n 

-=J>(r,(cJ))-r) X^AA(0,S), 



(6.19) 



Sij = —(rilr^+riTj- ±((ri\S(l) rj ) + ( rj \S(l) ri )) (6.20) 
is definite positive, where 5(1) denotes the reduced resolvent of P at 1. 
Remark 6.7 An alternative formulation for £ is 

S - = ^((r l |(lQ-PQ)- 1 (lQ + PQ)r,) + (r J |(lQ-P Q )- 1 (lQ + P Q )r i )), (6.21) 

where the projector Q and the operator Pq are defined in the spectral decomposition of P 
P = 2d\u )(u \ + P Q , with P Q = QPQ and Q = Q 2 = Q* = I - 2d\u )(u \. (6.22) 

Proof: Our assumptions imply that P is irreducible, doubly stochastic and that uq is 
invariant for P and P* . This together with (6.18) allows us to apply Thm 20.1 of [7] and 
the remarks p. 177 or the results of [25]. It remains to compute the covariance matrix. Let 
us define the centered random vector 

r(r(w)) = r(r(w)) - r (6.23) 

such that E Uo (f(ro)) = 0, where E„ denotes the expectation with invariant initial measure 
no- The first mentioned reference yields the following expression for the covariance matrix 

oo 

= E uo (r(r )if(To)j) + ^E uo (f (r );f (t^uJ)^ + r(T k (u))if(T )j), (6.24) 
k=i 
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for i, j = 1, 2, . . . , d, where f(r)j denotes the j th component of f(r) G 1, d . We compute for 
any k G N 

E^C^ro)^^),) = I^IPV,)-^, (6.25) 

Note that the right hand side of (6.25) is equal to 

^^{r i \P k (r j -2dr j u ))^ with (u \{r j -2dr j u )}=0. (6.26) 

By (6.22), for any v,w 6 

i> — 2dvuQ = and (iz;|Qv) = (u>|i>) — 2dwv. (6.27) 
Thanks to (6.26), we can write 



oo 



^P fc (r, - 2df jUo )) = J2 P Q^ ~ 2 dr jU0 )) = (Iq - Pq^Pq^ - 2dr jUo ), (6.28) 

k=l k=l 

where Iq is the identity reduced to the subspace QC 2d . Therefore, 

(I Q -P Q )-ip Q ^-(S(l) + I Q ), (6.29) 

where S(l) = QS(1)Q = (Pq — Ig^Q^d denotes the reduced resolvent of P at 1. Hence 

Sij = ^i\rj) - rirj - ^({ri\(S(l) + I Q )( rj - 2dr jUo )) 

+(r j \(S(l)+I Q )(r i -2dr i u ))) 

= -^(rilQrj) - Y d UnWrj) + <^(i)n» 

= ({n\{Q + 2S(l)) rj ) + (rj\(Q + 2S(l))n)) 

= ^((r i |(I Q -P Q )- 1 (IlQ + J PQ)r,) + (r J |(lQ-PQ)- 1 (lQ + PQ)r i )). (6.30) 



The convergence of Sj=i( r ( r i( u; )) ~~ r ) f° r an y initial measure po s.t. Po(t") = |or| 2 in 
distribution to ~ A/"(0, S) implies the convergence of the characteristic function and of 
its derivatives, which are continuous functions of the random variable 
In particular 

-d yj d y ^ iV ^n o (4=) |«=o = E l^,| 2 ^ fE(r(7l(5;)) -r),- f>(^)) - r) fc ) , 

(6.31) 

whose limit, as n — > oo yields the elements of the random diffusion matrix. 
We have 



29 



Corollary 6.8 Under the assumptions of Theorem 6.6, the following random variables 

converge in distribution as n — >■ oo: 

The random rescaled characteristic functions 



e -W^o (y/ ^) = la^pfe^^^ 1 ^^)^^^, (6.32) 



r ei± 

e iyX " , where X w ~ A/"(0, E), (6.33) 



one? t/ie random diffusion constants 



/ n n \ 

E i i 2 — Ew^)) - ^ E( r ^)) - T ^ — ► D ^ ( 6 - 34 ) 

() e/± \z=i i=i / 

E / ( T o\ R o(y/Vn - v, v)t ) dv- ( V(r(r;(^)) - r)j V(r(r ; (cJ)) - r) k ) 



(6.35) 



n,el± ■ 

where B^ fc is distributed according to the law of X^X^, where X u ~ A/"(0, E). 
Remark 6.9 in particular we get 

E W (B^) = E jfe . (6.36) 
Proof: For the case of initial density matrix po , it is enough to note that 

V f (T \R (-v,v)To)dv= [ TrRo(-v,v)dv = $o(0) = l. (6.37) 

One concludes using the convergence results stated in [7], p. 28 and 30. S 

At this point one may wonder if the assumption P = E M (C T (u;)) is enough to apply The- 
orem 3.10 and compare the results concerning the averaged distribution w(n) = E(W"(n)). 
The next proposition answers this question positively 

Proposition 6.10 Under the hypotheses of Theorem 6.6, assumption S' holds. Moreover, 
the diffusion constant D(t>) given by Theorem 3.10 is independent of v € T d . 

Proof: We need to consider 

M(y,y') = D(y,y')£ = ^ e i{yr{r)+y'r{r')) | T T /^ T r /| fj (g_ 38 ) 

r,r'e/ ± 



where £ = E^C^ <8> C w ), with f and C(w) defined above (6.11). We first observe that the 
M* (Y )-cyclic subspace generated by X, is given by X = span{|<r (g> a)} a ei±- Indeed, 
*i £l and X is invariant under D(Y) for any Y = (y, y 1 ) £ T M since D(y,y')\a a) = 
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e i{y+y')r{v)\ a ® a ). Then, for any C(vr,0) = E re / ± |vr(r))(r| and any {a a } aeI± , 

a a G C, we compute 

C(7r,0) C(tt,0) ^a CT |a(»cj) = ^ a CT |vr(cr) ® vr(o-)) 

= C(vr)0CV) ^2a a \a®a). (6.39) 

o-e/± 

This shows that I is invariant under (C w <8> C w )*, and thus under its expectation as well, 
which is enough to prove the claim. Moreover, (6.39) shows that, when restricted to I, any 
matrix C(n, 0) (8) C(7r, 0) acts like C(7r) (g) C(7r) does. Consequently, we can consider the 
finite measure \i on 6 2 d defined by (6.11) instead of the original measure v. Altogether we 
get 

M(-v,v)\ x = D(-v,v)\ x ^)C{k)®W)\t 

= »(K)C(ir)®C(ir)\i = £\ I , (6.40) 

which is independent of v G T d . It remains to show that *i is the only invariant vector 
under £\%. But the equation for the coefficients a a 

^2 M 71 ")^ 71 ") ® C(tt) ^ a a \a ® a) = ^ a a \a <g> a) (6.41) 

Tre&2d <r£l± a€l± 

is equivalent to X^7re6 2d M 71 ") ^- 1 ^) = a a, for all a. This is in turn equivalent to requiring 
that the vector a G C 2d with components {a cr } cr< =j ± be invariant under P T . P being 
irreducible by assumption, the only invariant vectors have constant components and thus 
are proportional to H 

Remark 6.11 The assumption S doesn't always hold for the case under study. In case 
= 0, the vector \l/ = ^(tt')67 2 I t ® r ') * s a ^ so invariant under A4(0, 0) = £ . 



We close this section by providing a general relationship between the diffusion matrix 
D w computed by means of W u (n) and the diffusion matrix D computed by means of the 
averaged distribution w(n) = E^W" '(ra)) in Theorem 3.10, provided both exist. 

To deal with the whole diffusion matrix at once, we let z G M d \{0} be fixed and consider 
the non-negative random variables 

vui _ ¥ 

D = (z\Bz), D% = (z\Y?)(Y?\z), where Y£ = n . (6.42) 

With these notations, K w ^(D n ) = E w Eww — > D by Theorem 3.10 and E W w( n )(D^) - 
D w where D w = (z|ID) w z) as n -)• oo, if the limit exists. 
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Proposition 6.12 Consider the initial condition = <po<gi\0) and assume the hypotheses 
of Theorem 3.10 hold for the distribution w(n). Further assume the random variables defined 
in (6.42) satisfy K W w^(D^) — > D u in distribution. Then 

E„(D") = lim EuEwu^DZ), (6.43) 

which implies E W (B W ) = D. 

Proof: Let = ¥,^uf n \(D^) > which converges in distribution to D u . By Theorem 5.4 
of [7], if sup n E w ((d^) 1+e ) < 00, for some e > 0, then the limit and expectation commute: 
lim^oo E w (d£) = E W (D"), which yields the result. Now, Remark 3.11 below Theorem 3.10 
implies the condition with e = 1. M 



Remark 6.13 When applied to the case discussed in the section, the proposition yields 

E W (W) = 'E= [ Bdv = H). (6.44) 

6.3 Large Deviations 

We can complete the picture for the model at hand by looking at its large deviations 
properties. Because the analysis reduces to the study of a finite state Markov chain, a large 
deviation principle is true for the model, see [13], Theorem 3.1.2. 

For A G M. d , let G M2d(K + ) be the matrix whose non- negative elements read 

n A ((7,r) =P(a,r)e <A|r(r)) , (6.45) 

where P is the transition matrix and r is the jump function. As P is irreducible, Tlx is 
irreducible as well. Hence, by Perron-Frobenius theorem, for all A G M. d , the largest real 
eigenvalue of Tl\, p(A) > 0, is simple and a (Tlx) C D(0,p(X)). As a function of A 6 M d , 
p(A) is real analytic. For every x G W 1 , define 

I(x) = sup «A|x> - ln(p(A))). (6.46) 

xeR d 

The function / is a good rate function. 

Theorem 6.14 Under the hypotheses of Theorem 6.6, the random variable Z n = ^ Y^=i r ( T j) 
satisfies a large deviation principle with convex good rate function I: For any a G I± and 
any T C R d 

- inf I(x) < lim inf — ln(P CT (Z n G T)) 

x <=r° n->oo n 

< limsup-ln(P CT (Z n G T)) < — inf I(x), (6.47) 
n xav 

where P CT refers to the initial law Pq(t) = 8 ajT 
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6.4 Example 

We end this section by a simple example which allows us to make explicit all quantities 
encountered so far. 



In dimension d = 1, we consider the jump function r defined by r(±l) = ±1, so that 
r = 0. We define the distribution /j, on ©2 by the Bernoulli process which assigns probability 

p > to the identity matrix and q = 1 — p to the matrix ( ^ J. The corresponding 



transition matrix reads 



P = ( P *) s.t. (P - z)"i = l*°><*°' + -M^L, (6.48) 
\9 JV l-z (p-q)-z 

with * = ^ (T), Y>i = 75 f \Y Consequently 5(1) = g^j-. Thus, the averaged 
diffusion constant B = £ computed from (6.20) with r = ^ reads 

X = -\{r\r)-±2(r\S(l)r) = Z. (6.49) 

Therefore the random variable %"(n)(^) converges in distribution to X w ~ M(0,p/q) 
and the corresponding random diffusion constant B w is distributed according to the chi- 
square law with density |/(-|) where 

e -t/2 

V 27rt 



Next we compute the rate function I such that P (7 (E^(n)( „ ) G T) ~ e -« inf -er /(x) ) for 
large. With our choice of jump function, the matrix II A reads 



n 



Ux=( Pe x ) s.t. detn A =p-q and Tr n A = 2pcosh(A). (6.51) 



qe pe 



Thus we have 



p(X) = pcosh(A) + yp 2 sinh 2 (A) + q 2 (6.52) 
so that sup AeR (xA — ln(p(A))) is reached at X(x) = arsinh ^ — /f=f ) f° r M < 1- Hence 



qx \ ( q + y 7 ^ 2 + ^ 2 (<? - p) 



p\/l — x 2 J \ \/l — x 2 



I(x) = x arsinh j - In j ' v ' ^ " ^- j if |.r| < ! ((>.o.:i) 

and I(x) = oo otherwise. 



33 



7 Generalization 



In the light of the last example, it is natural to generalize the results of Section 3 to the 
case where the random coin matrices C(cu) are distributed according to a Markov process, 
in the spirit of [29], [17], [15]. We briefly do so in this last section, mentioning the main 
modifications only and considering finitely many coin matrices for simplicity. 

Consider a finite set {Ci, C2, ■ ■ ■ ,Cf} of unitary coin matrices on C 2d and assume that 
for any n G N, the set of random matrices {C(u)\), • • • , C(ui n )} is determined by a Markov 
chain characterized by an initial distribution {po{j)}f=i and an irreducible transition matrix 
{P(j, k)}j^ e u... Correspondingly, for any Y G T d x T d , the sequence of matrices 
{Muj (X)}jeN has the same distribution, assuming the Cj <g> C k are distinct, so that 

P({M Wn (y), M Un _, (¥),■■■ M m (Y)} = {M kn (Y), M^Y), ■ ■ ■ M kl (Y)}) 

= po(fci)P(fci, k 2 ) ■ ■ ■ P(A; n _i, k n ). (7.1) 

We introduce the matrices 

M jk (Y) = P T (j, k)M k (Y), j, k G {1, ■ ■ ■ , F}, (7.2) 
acting on C 2d (g C 2d and the operator acting on C F <g (C 2d <g C 2d ) 

M(Y)= Yl \ej)(e k \®M jk (Y), (7.3) 
j,ke{i,-,F} 

where the e/s are the canonical basis vectors of C F and the "bra"s and "kef's refer to the 
usual scalar product in C F . Similarly, let Mq(Y) be the vector in M 4d 2(C) F defined by 

Mo(Y)= Yl Poti)\ej) ® Mj(Y). (7.4) 

jG{l,-,F} 

If xi = J2je{i ■■■ F} \ e j) G ^ F anc ^ " i s * ne identity operator on C 2d (g C 2d , we obtain 

EiM^^M^iY) ■ ■ ■ M m (Y)) = ( Xl 0l|M(yr- 1 M o (y)). (7.5) 
Hence, with = ^ re / ± \t r) and <&o £ < C 2d ® C 2d , we can write 

E(COO) = (Xi ® *i|M(y)"- 1 M (y)$ ), (7.6) 

where xi ® *i and M (y)$ = EjPoC?)l e j) (g Mf(y)$ belong to C F <g) (C 2d C M ) and 
"bra"s and "kefs should be interpreted accordingly. 

This brings us back to the study of large powers of an operator, M(y), which has 
essentially the same properties as M(Y): 

Lemma 7.1 The matrix M(Y) acting on C F (g (C 2d <g C 2d ) is analytic inY G C 2d x C 2d . 
Assume P is irreducible, and let Xp G C F be the unique real valued vector s.t. P T Xp = Xp 
and (Xplxi) = 1- Then for any v G T d , 

M{-v,v)\x P ® *i) = |Xp® *i) and M*(-7;,t;)|xi (8) *i) = |xi ® *i>- (7-7) 
For any y G T d x T d ; 

SprM(y) = SprM*(y) < 1. (7.8) 
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Proof: The first two identities follow from explicit computations. The second prop- 
erty is a consequence of the fact that if one endows <8> (C 2d £3 C 2d ) with the norm 
|| ej <£> VPjHoo : = maxj \\ C 2d^ C 2d, then M*(Y) becomes a contraction. This is due to 
the fact that P, as a stochastic matrix, is a contraction with the sup norm on C F . The 
spectral radius of M* (Y) thus cannot exceed one. H 

We shall work under the 

Assumption S": For all v G T d , 

a(M(-v,v)\i*) n dD(0, 1) = {1} and the eigenvalue 1 is simple, (7.9) 
where X* is the M* (Y)-cyclic subspace generated by %i ® 

Then we can proceed with a spectral analysis similar to that of Section 3. Ignoring the 
restriction |j* in the notation, we note that Assumption S" implies that for all v 6 T d , a 
complex neighborhood of T d , we can write 

(M(-v,v) - z)- 1 = JL- + S v {z), z g a(M(-v,v)), (7.10) 

where P = \\Xp ® ® *i| is independent of v, and the reduced resolvent S v (z) has 

the same analyticity properties as that of M(—v,v). Moreover, introducing 

A(Y) = Y,\e J )(e J \®D(Y) (7.11) 

3 

we see that M(y — v, v ) = A(y, 0)M(— v , v) where 

A(y, 0) = J2 \ej)(ej\ ® (I + F 1 (y) + F 2 (y) + 0(\\yf)), (7.12) 

3 

see (3.24). Therefore, if (y, v) 6 ^(0, yo)xT d , for y > and v > small enough, Lemma 3.6 
holds for M(y — v , v). Applying the same perturbation formulas for the isolated eigenvalue 
of M(y — v,v), noted Xi(y, v) again, for y € C d small enough, we reach the same conclusions 
by explicit computations: 

Mv,v) = l + ^^yr^-^ylD^ + O^Hyll 3 ), (7.13) 

rel± 

for all v G T d . The first order term in y is the same as the one of (3.27). On the other 
hand, the explicit form of the quadratic term in y which defines the (analytic) diffusion 
matrix ©(f), depends on the transition matrix P 



(y\B(v)y) = - 
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Moreover, the corresponding rank one projector P{y, v) is analytic and thus tends to P 
uniformly in v G , as y — > 0. With Q(y,v) = I — P(y,v), the matrix Q(y,v)M.(y — 
v,v)Q(y,v) has spectral radius strictly small than one, for v G T„ and y G C d small 
enough. 

Therefore we can state that all conclusions drawn in Section 3, e.g. Theorem 3.10, and in 
Section 5, e.g. Theorem 5.2, for i.i.d coin matrices under Assumption S are true for finitely 
many coin matrices forming a Markov chain, under Assumption S" , mutatis, mutandis. 
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